GREEK LETTERS IN RANDOM STAIRCASE TABLEAUX* 
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Abstract. In this paper we study a relatively new combinatorial object called 
staircase tableaux. Staircase tableaux were introduced by Corteel and Williams 
in the connection with Asymmetric Exclusion Process and has since found in- 
teresting connections with Askey-Wilson polynomials. We develop a proba- 
bilistic approach that allows us to analyze several parameters of a randomly 
chosen staircase tableau of a given size. In particular, we obtain limiting distri- 
butions for statistics associated with appearances of Greek letters in staircase 
tableaux. 



Key words and phrases: asymmetric exclusion process, asymptotic normality, stair- 
case tableaux. 



1. INTRODUCTION 

An interesting combinatorial structure, called staircase tableaux, was introduced 
in recent work of Corteel and Williams [11, 12]. Staircase tableaux are related to the 
asymmetric exclusion process on an one-dimensional lattice with open boundaries, 
the ASEP. This is an important and heavily studied particle model in statistical 
mechanics (we refer to [12] for some background information on several versions of 
that model and their applications and connections to other branches of science). 
The study of the generating function of the staircase tableau has given a combina- 
torial formula for the steady state probability of the ASEP. Explicit expressions for 
the steady state probabilities were first given in [15]. In their work [12, 11] Corteel 
and Williams used staircase tableaux to give a combinatorial formula for the mo- 
ments of the (weight function of the) Askey-Wilson polynomials; for a follow-up 
work see [7]. 



f LIAFA, Universite Paris Diderot-Paris 7, F-75205 Paris, France, dasse@liafa.jussicu.fr. 

X Institute of Mathematics, Polish Academy of Sciences, ul. Sniadcckich 8, 00-956 Warszawa, 
Poland and Department of Mathematics and Information Science, Warsaw University of Technol- 
ogy, PI. Politechniki 1, 00-661 Warszawa, Poland, phitczenko@math.drexel.edu. 

* The work of the first author was carried out while she held an ANR Gamma internship at 
LIPN, Universite Paris Nord, under the direction of Frcdcrique Bassino (LIPN) and Sylvie Corteel 
(LIAFA). She would like to thank both of them for their guidance, the members of LIPN for their 
hospitality, and ANR Gamma for the support. The second author was partially supported by 
NSA Grant #H98230-09-l-0062. Most of his work was done during his stay at LIPN in July 2010 
and completed while he was at the Institute of Mathematics of the Polish Academy of Sciences 
and the Technical University of Warsaw in the Fall of 2010. He would like to thank Frcdcrique 
Bassino for the invitation and acknowledge the hospitality of these institutions. 



2 



SANDRINE DASSE-HARTAUT AND PAWEL HITCZENKO 



The authors of [12] called for further investigation of the staircase tableaux be- 
cause of their combinatorial interest and their potential connection to geometry. 
In this paper we take up that issue and study some basic properties of staircase 
tableaux. More precisely, we analyze the distribution of various parameters associ- 
ated with appearances of Greek letters a, j3, S, and 7 in randomly chosen staircase 
tableau of size n (see the next section or, e.g. [12, Section 2], for the definitions 
and the meaning of these symbols). 

Staircase tableaux are generalizations of permutation tableaux (see e.g. [6, 9, 
10, 20] and references therein for more information on these objects and their con- 
nection to a version of ASEP referred to as the partially asymmetric exclusion 
process; PASEP). For permutation tableaux, the authors of [6] developed a proba- 
bilistic approach that later allowed the derivation of the limiting (and even exact) 
distributions of various parameters of the permutation tableaux. Our goal here 
is the same: in Section 3 we develop a probabilistic approach parallel to that of 
[6] that allows us to compute generating functions of various quantities associated 
with staircase tableaux (see Corollary 6 in Section 4 and Proposition 11 in Subsec- 
tion 6.1.3). As a consequence, we obtain the exact or limiting distributions of the 
parameters we study. 

Our main results are gathered in Section 5 and may be summarized by stating 
that the five parameters we consider have asymptotically normal distribution when 
normalized in a usual way (i.e. centered by the mean and scaled by the square root 
of the variance). Thus, for example, Theorem 9 asserts that the number A n of a 
or 7 on the diagonal of a randomly chosen staircase tableau of size n has expected 
value n/2, variance (n+l)/12 and that (A n — n/2)/ y/n/12 converges in distribution 
to the standard normal random variable. We refer to Theorems 8 and 9 for precise 
statements. 

In Section 6 below (see comments at the beginning of Subsection 6.1.3 and a 
remark at the end of it) we find that one of the parameters we study coincides 
with a generalization of Eulerian numbers (see [23, sequence A060187]) related to 
Whitney numbers of Dowling lattices (see [23, sequences A145901, A039775] and 
[16, 2, 3, 5] for definitions and further information on these numbers). This rather 
unexpected and intriguing connection has not been explained and merits, perhaps, 
further studies. One consequence of our work is that the triangle of numbers [23, 
sequence A060187], when suitably normalized, satisfies the central limit theorem. 
As far as we can tell this result is new (although it is an easy consequence of a 
general theorem of Bender [1]). Limit theorems for a related sequence [23, A145901] 
are established in [5]. This link to the Whitney numbers of Dowling lattices may 
have unraveled connection to geometry alluded to in [12] as the sequences A060187, 
A145901, and A039775 from [23] all have a very strong geometrical flavor. 

2. DEFINITIONS AND NOTATION 

We recall the following concept first introduced in [11, 12]: A staircase tableau 
of size n is a Young diagram of shape (n,n — 1, . . . ,2, 1) whose boxes are filled 
according to the following rules: 

• each box is either empty or contains one of the letters a, (3, 6, or 7; 

• no box on the diagonal is empty; 

• all boxes in the same row and to the left of a (3 or a S are empty; 

• all boxes in the same column and above an a or a 7 are empty. 



GREEK LETTERS IN STAIRCASE TABLEAUX 



An example of a staircase tableau is given in Fig. 1(a). 
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Figure 1. (a) A staircase tableau of size is 7. Its top row is 
indexed by /3, the next one by a. (b) Its type (o o • • • o o). (c) 
The same tableau with it's and q's. 



As we mentioned, staircase (or earlier permutation) tableaux were studied in 
the connection with ASEP. Because of the importance of this connection we briefly 
recall its nature. The ASEP is a Markov chain on words of size n on an alphabet 
A = {o, •} consisting of two letters. Each such word represents an one-dimensional 
lattice of length n with some sites occupied by particles (represented by •), and 
others not (represented by o). A particle can only hop to the right or the left 
(with the probabilities u and q, respectively), provided that the adjacent site is 
unoccupied, or enter or quit the lattice. Entering from the left (right) happens with 
the probability a (resp., S) if the first (last) site is unoccupied. Exiting to the left 
(right) happens with the probability 7 (resp., (3) if the first (last) site is occupied. 
At a given time one of the n + 1 possible locations for a move is selected (uniformly 
at random) and, if possible, a transition described above is performed with the 
given probability. We refer to [14, 15, 17] or [12] for more detailed description and 
further references. 

To describe the connection to staircase tableaux, define the type of a staircase 
tableau S of size n to be a word of the same size on the alphabet {o, •} obtained by 
reading the diagonal boxes from northeast (NE) to southwest (SW) and writing • 
for each a or S, and o for each /3 or 7. (Thus a type of a tableau is a possible state 
for the ASEP.) Fig. 1(b) shows a tableau with its type. We also need a weight of a 
tableau S. To compute it, we first label the empty boxes of S with u's and g's as 
follows: first, we fill all the boxes to the left of a j3 with it's, and all the boxes to 
the left of a 6 with g's. Then, we fill the boxes above an a or a 5 with u's, and the 
boxes above a j3 or a 7 with q's. When the tableau is filled, its weight, wt(S), is a 
monomial of degree n(n + l)/2 in a, /3, 7, 5, u and q, which is the product of labels 
of the boxes of S. Fig. 1(c) shows a tableau filled with u's and q's. Its weight is 
a 3 /? 2 <5 3 7 3 uV. 

Corteel and Williams [12, 11] have shown that the steady state probability that 
the ASEP is in state a is 

Z a (a,f3,j,6,q,u) 



where Z n {a, (3, 7, 5, q, u) 



S of size n 



Zn{a,l3,j,5,q,u)' 

wt(S) and Z a (a, ft, 7, 5, q, u) 



S of type a 
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We denote the set of all staircase tableaux of size n by S ni n > 1. It is known 
that the cardinality of S n is 4™n!. There are several proofs of this fact (c.f. [7] for 
one of them and for references to further proofs). All these proofs are based on 
combinatorial approaches and we wish to mention that a probabilistic technique 
that we develop in this paper provides a yet another proof of that fact. We present 
it in Section 4 below as an illustration of how our method works. 

We now define some parameters that are the objects of our study. Let * be a 
subset of the set of symbols {a, (3, 5, 7}. We say that a row of a staircase tableau 
is indexed by * if its leftmost entry is a member of *. For the sake of brevity we 
will refer to rows indexed by * simply as * rows. Thus, for example, the number of 
a/7 rows is the number of rows indexed by a or 7. The tableau in Figure 1(a) has 
two a/7 rows, the second from the top (indexed by a) and the bottom (indexed by 
7). For a given staircase tableau S <E S n we denote this quantity by r n (S) and we 
occasionally skip the subscript n if there is no risk of confusion. As we demonstrate 
in Section 3 below this parameter plays a fundamental role in our approach. Other 
parameters we consider are: the total number of entries (3 or 5 ((3/5 for short), the 
total number of entries a or 7 (a/7), the number of entries (3/5 on the diagonal of 
the tableau, and the number of entries a/7 on the diagonal. For a given tableau 
S e S n these parameters will be denoted by A n (S), T n (S), B n (S), and A n (S), 
respectively. We gather the names and the notation we use for the parameters we 
consider along with their values for a tableau given in Figure 1(a) in the following 
table 



parameter 


notation 


value in Fig. 1(a) 
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An 


5 


total # of a/7 


r 


6 


# of (3/5 on diagonal 




3 


# of a/7 on diagonal 
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4 



As we mentioned earlier our viewpoint is probabilistic. Thus, we equip the set 
S n with the uniform probability measure denoted by P„. This means that for each 
S € S n we have 

p ^ = ^: . 

As is customary we refer to a tableau chosen according to that measure as a random 
tableau of size n. We denote the integration with respect to the measure P„ by 
E„. Our goal is to analyze probabilistic properties, like expected values, variances, 
and exact or limiting distributions, of random variables (or statistics) r n , A„, T n , 
B n , and A n . We follow the probability theory conventions of [22], and the reader 
is pointed there for any unexplained terms. 

3. Preliminaries and outline of the argument 

In this section we detail main ideas beyond our approach and also derive basic 
properties of the fundamental parameter, i.e. the number of a/7 rows. 

Our method is analogous to what has been done in the case of permutation 
tableaux (see [6] or [20]). Let us recall at this point that permutation tableaux 
have been used to give a combinatorial description of a stationary distribution 
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for the PASEP. We refer to e.g. [9, 10, 6, 20] for the definition, connections to 
PASEP, further properties and details. Just as PASEP is a particular case of 
ASEP, permutation tableaux of size n are in bijection with a subset of staircase 
tableaux of that size corresponding to the case 7 = 5 = 0. 

The approach used in [6, 20] for the permutation tableaux was to identify a 
fundamental parameter, trace its evolution as the size of a tableau is increased 
by 1, and then use successively conditioning to reduce the size of a tableau. Wc 
refer the reader to either [6] or [20] for more details, here we only recall that this 
fundamental parameter was the number of unrestricted rows, U n , in a permutation 
tableau, and that its conditional distribution was given by 1 + Bin(J7„_i, 1/2) (here 
Bin(n,p) denotes a binomial random variable with parameters n and p). This is 
to mean that if a size of a permutation tableau with U n _i unrestricted rows was 
increased from n — 1 to n, then the number of unrestricted rows in this extension 
had (the conditional) distribution 1 + Bin([/„_i, 1/2). As it turns out, in the case 
of staircase tableaux the role of a fundamental parameter is played by the number 
of a/7 rows. 

To make our approach work we need to know two things. One is the knowledge of 
the conditional distribution of statistics of interest in extension of a given tableau of 
size n—1. This would allow us to use the so-called tower property of the conditional 
expectation to reduce the size of staircase tableaux by 1. Roughly speaking, this 
reducing would amount to relating a function of a statistic on, say, S n to another 
function of the same statistics on iS„_i. Keeping track of how the functions are 
related would in principle allow for an iteration of this process to reduce the size 
of the tableaux to 1. There is, however, another difficulty. Namely, when passing 
from S n to iS„_i, not only a function of statistic changes but also a measure over 
which we integrate. More specifically, the image of the uniform probability measure 
P„ on S n under a natural mapping associated with our reduction does not yield 
the uniform probability measure P„_i on <S„_i. To handle this difficulty wc need a 
second ingredient which we refer to as the change of measure. It simply describes the 
relation between the two probability measures and allows to return to the uniform 
probability measure on S n -i after the reduction (see Section 3.2 for more details). 

In the reminder of this section we develop the details for the number of a/7 rows, 
a key statistic in further considerations. We do it in two separate subsections. In 
the first we analyze the evolution of the number of a/7 rows as the size of the 
tableaux increases. The key step is the computation of the conditional generating 
function of r„ as the size of the tableau is extended from n — 1 to n (Lemma 3 
below). In the second subsection we discuss the change of measure. 

3.1. The conditional generating function for r n . Our analysis progresses in 
steps. We first find the number of extensions of a given tableau of size n — 1. Wc 
then find the conditional distribution of the number of a/7 rows. This enables us 
to compute the conditional generating function of the number of such rows in an 
extension of a given tableau of size n—1. The formula expresses this generating 
function in terms of the number of a/7 rows of the tableau being extended and is 
of crucial importance for our approach. 

To begin the analysis we need to briefly recall the evolution process of staircase 
tableaux described by Corteel, Stanton, and Williams in [8]. Let S be a staircase 
tableau of size n—1. To extend its size by 1, we add a new column of length n at 
the left end and we fill it according to the rules. Any such filling gives a tableau 
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of size n which is an extension of S. The number of extensions clearly depends on 
the number of a/7 rows of S and is as follows: 

Lemma 1. Let S £ S n -\ be a tableau with r n -\ = r n -i(S) a/7 rows. Then there 
are 4 • 3 r,l_1 different staircase tableaux of size n that are extensions of S . 

Proof. We consider two cases: 

1. If r„_! = then all n — 1 rows of S are (3/8 rows and hence the top n— 1 boxes 
of the new column have to be empty since no entries are allowed to the left of a 
(3/8 in the same row. Thus, if r„_i = we obtain four different tableaux of size n 
by putting one of the four symbols in the bottom box of the new column. 
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Figure 2. Extension of a tableau of size 6: (a) a/7 is put in the 
bottom corner of the 7th column. All boxes above it are to be 
empty, (b) (3/8 is put in the bottom box and a/7 in the 3rd of the 
four a/7 rows in the 7th column. A box in an a/7 row above is to 
be empty; |fs in two a/7 rows below it indicate that these boxes 
can be filled with (3/ 8 or left empty. 



2. If r n -\ > 1 we count the number of extensions as follows (see Figure 2 for an 
illustration): we can either put one of the symbols a/7 in the bottom corner (and 
then we are forced to leave all other boxes in that column empty), or we can put 
a, (3/8 in the bottom box of the new column. In the latter case, we need to fill the 
r„_i boxes in the new column corresponding to the a/7 rows of S. According to 
the rules, if we put an a/7 in any of them, then we need to leave all boxes above it 
empty, otherwise we have a complete freedom. It follows from these considerations 
that any tableau of size n with r n _i a/7 rows gives 4 • 3 r " _1 different staircase 
tableaux of size n. Indeed, if r„_i = then there are 4 extensions and if r„_i > 1 
then there are 

2 + 2(2(1 + 3 + • • • + 3 r "- 1_1 ) + S''— 1 ) 

extensions. Here, the first 2 is from putting an a/7 in the bottom box of the new 
column, the next 2 is from putting a /3/8 in that box, the term 2-3 4 " 1 , 1 < i < r n _i, 
is from putting the first a/7 in a box of the new column corresponding to the zth 
(counting from bottom) a/7 row of S as then there are 3 l_1 ways of filling the 
earlier i — 1 boxes with symbols (3, 8, or leaving them empty, and finally, the 3 r " _1 
term comes from not putting an a/7 in any of the r n _i boxes and thus filling them 
with /3/S's or leaving them empty. 
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Summing the above gives 

2 + 2 ( V'"^ - 1 + 3 r -^ = 2 + 2(3 r "" 1 - 1 + S''— 1 ) = 4 • 3 r »-\ 

as claimed. □ 

As the next step we find the conditional distribution of the number of a/7 rows 
over the extensions of a given tableau of size n—1. To do that we phrase the preced- 
ing discussion in a more probabilistic language using what Shiryaev [22, Chapter I, 
§8 ] refers to as "decompositions" (which is just a special case of conditioning with 
respect to cr-algebra). 

Note that every tableau from S n is an extension of a unique tableau S from <S„_i. 
Therefore, denoting by D$ the set of all tableaux from S n which are obtained from 
S by the process described above (as we just discussed, there are 4 • 3 r «-i( s ) such 
tableaux), we can write 

S n = (J D S , 

where (£ ) s)se5„_i are pairwise disjoint, non-empty subsets of S n . We denote this 
decomposition of S n by T> n -\. We wish to compute the conditional probabilities 
P( • \D n _\) and the conditional expectations E( • \D n _i) with respect to this de- 
composition. We have 

Lemma 2. Let S € S n -\ be a tableau with r„_i = r n -\{S) a/7 rows and let r be 
the number of such rows in any of its extensions to a tableau of size n. Then 

P(r = r n _ 1 + l\D s ) = Y -± F — i , 

0, l,...,r„_i 

Proof. Let S and r be as in the statement. Clearly, the possible values for r are 
r n -i + 1, r„_i, r n -i — 1, . . . , 1, and we need to know the probability of each of 
these possibilities. First, r = r„_i + 1 means that we must have placed an a/7 
in the bottom box of the new column since this is the only way we can increase 
the number of a/7 rows. Consequently, all other boxes in the new column are to 
remain empty. Obviously there are two such extensions which in view of Lemma 2 
gives 

P(r = r„_! + 1\D S ) = ^-J-t- = ^-7, 

which justifies (1). 

Next, for k — 0,1,... , r„_i we compute P(r = r n _i — k\Ds). Since k is the 
number of /3/S's that we put in the r„_i allowable (i.e. corresponding to a/7 rows 
of S) boxes, r = r„_i — k means that we put a (3/5 in the bottom box and additional 
k /3/<5's in the r„_i allowable boxes above it. If we do not put an a/7 in any of 
those k boxes, we have 2 • 2 fe ( r " fc _1 ) possibilities (2 is for putting either a [3 or a S 
at the bottom, and the rest accounts for putting /3/S's in any k of the allowable 
r n _i boxes). If we do put an a/7 in one of the boxes we need to pick k + 1 of 
the allowable r„_i boxes, put an a/7 in the topmost of them and put /3/5's in 
the remaining k of them (and at the bottom). This gives 2 fe+2 ('>_r 1 1 ) possibilities. 



(1) 

and, for k 
(2) 
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Adding up the two pieces and dividing by the total number of extensions given in 
Lemma 2 we obtain, 

p , kW) - 2 ' 2 *CV) + 2 fc+2 to) _ 2 fc+1 (2fe) + (V)) 

as required. □ 

Lemma 2 completely describes the conditional distribution of r n given T) n _\ and 
leads, in particular, to the following basic relation: 

Lemma 3. For a complex number z and n > 1 (with the understanding that r = 0) 

(3) E(* r "|Z>n-i) 1 : ' 

Proof. Recall that as we extend a tableau with r n _i a/7 rows, we have r„ = r„_i+l 
if and only if we put an a/7 at the bottom of the new column and r n = r n _i — k if 
and only if we put a /3/S there and additional k /3 / <5's in the r n _i boxes above it. 
Therefore, letting C* denote the event that we put the symbol * at the bottom of 
the new column we have 

E(* r »|E>„_i) = E(z r «I Ca ^\T> n -i) + E(z r «I Cli , t \'D n -i) 

= ^ n " 1+1 ^T7 + E 2 r "- 1_fc P(rn = r„_i - fc|2>„_i) 



2 13/ 4-3^-1 ^ V. V./.' • 1/ k 



r„_i-l 



The last term is 

r„_i-l 



- — y 2 fe +v»- fe - 1 r*- 1 ) = — ^ — ((2+*) r »- 1 -^"-^ 



2 . 3r„_i Z^ v k + 1 I 2 • 3 r "-i 

k=a v 7 

so that we obtain 

E(z r " |2? n -i) = 2 ( 3 ) +2 (^j +2 (^j - 2 ( 3 ) 



l + z /2 + z x rn " 1 



This proves (3). □ 

As we already mentioned, Lemma 3 (or its versions) plays a crucial role in our 
approach. 

3.2. The change of measure. In this subsection we describe the second ingredi- 
ent, namely the change of measure. It it is necessitated by the fact that there are 
two different probability measures on <S„_i that naturally appear in our consider- 
ations. The first is, of course, the uniform measure P„_i. We discuss the second 
one and the relation between them below. 

Consider <S„_i, the set of all staircase tableaux of size n — 1. The second prob- 
ability measure on <S„_i is obtained from the uniform probability measure on S n 
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by "collapsing" all the elements of S n that are extensions of the same element 
S e <S n -i- More formally, consider a map / : S n — > 5„_i defined by /(T) = 5" if 
and only if T is an extension of S. The measure of interest is the image of P„ under 
/. We denote this measure P n (S) (there is an apparent ambiguity of notation here, 
however, it disappears once we remember whether S is in S n or <S„_i). Since both 
of these measures appear in the course of our argument, it is important to find the 
relationship between them. But this is straightforward: since a tableau S G <S„_i 
with r„_i a/7 rows gives 4 • 3 r ™- 1 tableaux in S n we have 

(4) P ( S ) 4 ' 3r "' 1 - 4 - 3r '" 1 l^-il 1 _ 4 r „_J^-il p 

(4) P " (6) ~ |S„| - \S n \ \S n ^\ " 4 3 \S n \ Pn - l{S} - 

Consequently, for any random variable X n _i on <S„_i we have 

(5) E„X„_! = E n ^4-y"- ll -^f^X n ^ = 4 %^E„_ 1 3'»- 1 X n _ 1 . 

|<->n| |<->n| 

Here we have used the same convention as above; for a random variable X on <S„_i, 
E„_iX denotes the expectation with respect to the uniform measure on while 
E„X denotes the expectation with respect to the measure that is induced on 5„_i 
by the uniform measure on S n . 

The relations (3) and (5) are key and will allow us to analyze the distributions 
of the various statistics on S n . Note that (4) and (5) are true regardless of whether 
we know the cardinalities of S n -i and S n or not. As a matter of fact, one can use 
(5) to provide a yet another argument that |<S„| = 4™n! as we will see in the next 
section. 



4. Generating function for the number of a/7 rows and some 

CONSEQUENCES 

In this section we demonstrate how we intend to apply (3) and (5) to derive 
recurrences for generating functions that upon solving yield information on the 
corresponding statistics. We focus here on the number of a/7 rows which, on one 
hand is the easiest to analyse, and on the other is central in the analysis of other 
statistics. 

Proposition 4. For every complex number z we have 

on n 

(6) E n z r " = — l[(z + 2k-l). 

^"I fc=i 

Proof. By the basic properties of the conditional expectation (see e.g. [22, For- 
mula (16), p. 79]) the expectation on the right-hand side is equal to E n E(z Trl |X>„_i). 
Using first (3) and then (5) we further get 

Ez r n _ Z + l (Z + 2V-- 1 _ (s + l HISn-J 



\$n\ 



2{z + l)^^E n ^{z + 2Y-K 

<->n 



Applying the same procedure with z replaced by z + 2 and n by n — 1 we obtain 



E n z r - = 2 2 (z + l)(z + 3)%r£n- 2 (^ + 4) 



r n -2 
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Further iteration yields 

E„z r " = 2 n ~\z + l)(z + 3) . . . (z + 2n - 3)^\ei(z + 2(n - l)) ri . 

Since |<Si| = 4 and Ei(z + 2(ra— l)) ri = |(z + 2(n — 1)) + | the above can be written 

as 

2 n 

E„z r " = — (z + l)(z + 3) . . . (z + 2n - 1), 

|«->n| 

which is precisely (6). □ 

Proposition 4 has a number of consequences. First, by putting z = 1 in (6) we 
obtain an independent confirmation of the count of staircase tableau of a given size. 

Corollary 5. Let S n be the set of all staircase tableaux of size n>l. Then 

\S n \=4 n n\ 

Note that once this corollary is known (4) and (5) simplify to 

(7) P n (5) = — P„_i(S) and E„X„_i = -E„_i3 rn_1 X„_i, 

n n 

respectively, and this is the form we will be using from now on. 
Next, by combining this corollary with (6) we obtain 

Corollary 6. The probability generating function of the number o/a/7 rows in a 
random staircase tableau of size n is given by 

k=l fe=l 

The last corollary gives, in turn, a complete information on the distribution of 

r n . 

Corollary 7. For every n > 1 we have 

n 

(8) r„ = ^J fe , 

fe=i 

where Jk 's are independent and Jk is a random variable which is 1 with probability 
l/(2fc) and with the remaining probability. In particular, 

( 9 ) E ^ = Ei = ^< Var(r " ) = E 2l( 1 -2l)=^-^' 

fc=i fc=i v 7 

where H n = Y^ik=i \ an< ^ ^ = Sfe=i W are harmonic numbers of the first and 
second order, respectively. Furthermore, 



_ Inn 

(10) "2 JV(0, 1). 

/ 1 R n 

V 2 

Proof. Note that a factor 



z + 2k-l =± + 

2fc 2fc 1 2fc j ' 

given in the previous corollary is the probability generating function of a random 
variable Jk which is 1 with probability 1 / (2k) and with the remaining probability. 
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Since the product of the probability generating functions corresponds to adding 
independent random variables, we obtain (8) and thus also (9). Finally, since 
Jfc are uniformly bounded and variances of their partial sums go to infinity, the 
Lindeberg condition for the central limit theorem (see e.g. [22, Chapter III §4]) 
holds trivially. Since as n — » oo, E„r„ ~ var(r„) ~ ^j 1 , (10) holds. □ 

5. Main results 

Our technique allows us to obtain further results concerning the distributions 
(sometimes exact, sometimes only asymptotic) of the statistics discussed above. 
We gather our results in the following two statements, concerning the total number 
of entries and the number of entries on the diagonal, respectively. Recall that A„ 
and T n denote the total number of f3/5 and a/7 in the tableau of size n, respectively. 

Theorem 8. Consider the set S n with the uniform probability measure P n . Then: 
(i) For every n > 1 we have 



(11) A n 'Vt 



fc=i 

where (Ik) are independent and P(Ik — 1) = 1 — j^, P(Ik = 0) = ^ . In 
particular, 

TT TT ff(^) 

(12) E n A n = n - — , var(A„) = — — , 



and, as n — > oo, 



A„ - n + \ In n 



(13) 2 ^JV(0,1). 

y'llnn 

(ii) For every n > 1 

T = A 

In particular, (11), (12), and (13) hold with A„ replaced by T n . 

Our second result concerns the entries on the diagonal. Recall that A n (resp. 
B n ) denote the number of a/7 (resp. (3/S) on the diagonal of a tableau of size n. 
For these parameters we get 

Theorem 9. The expected value and the variance of the number A n ofa/j on the 
diagonal of a random staircase tableau of size n are, respectively, 

TL Tl ~\~ 1 

(14) E n A n = - and var(A„) = 
Moreover, 

(15) ^-^AiV(0,l). 

Furthermore, for every n > 1 we have 

(16) B n I A n . 

In particular, (14) and (15) hold for B n in place of A n . 
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Remark: While it seems intuitively clear that the expected number of letters a/7 
on the diagonal is n/2 as (14) asserts, the expression for the variance is much less 
intuitive. It implies, in particular, that the entries a/7 and ft/ 8 along the diagonal 
are not chosen independently from one another with equal probabilities as one 
might have hoped (if that were the case the variance would be n/4). 

6. PROOFS 

We begin by observing that, as is obvious from the definition, staircase tableaux 
are symmetric under taking the transpose and exchanging the a's with the /3's and 
the 7's with the (5's. It therefore is immediate that parts (ii) of both theorems 
follow from the respective parts (i). Furthermore, the proof of Theorem 8 (i) may 
be completed by using Corollary 7 and the relation 

(17) r n (S) + A n (S) = n 

since then 

n n 

A n = n - r n = ^(1 - Jfc) = ^ I k - 

fe=l k=l 

But (17) is clear once we notice that for any row of a staircase tableau exactly one 
of the following statements is true 

• it contains a 0/5 

• it is indexed by a/7. 

This proves Theorem 8 and we now turn our attention to a considerably more 
involved proof of Theorem 9 (i). 

6.1. Proof of Theorem 9 (i). The idea is the same as for Corollary 7 except that 
this time we will actually need the joint probability generating function of A n and 
r n . The final expression will turn out to be substantially more complicated than 
what we encountered earlier and thus harder to analyze. Nonetheless, the situation 
is quite analogous to the case of the number of rows in permutation tableaux (see 
[20, Section 4]). For the reader's convenience we break up our proof into several 
steps each of them discussed in a separate section below. We now briefly outline 
the major steps in the proof indicating the section in which they are treated. We 
begin in the forthcoming section with the derivation of the probability generating 
function. Its coefficients satisfy certain recurrences. This, in particular, enables 
us to derive the exact formulas for the expected value and the variance of A n 
(see Subsection 6.1.2 below). Furthermore, the nature of these recurrences suggests 
that the coefficients are related to the classical Eulerian number associated with the 
number of rises in random permutations. In fact, our coefficients exactly match the 
numbers often called the "Eulerian numbers of type B". We establish and discuss 
further this connection in Subsection 6.1.3. We think it is of independent interest 
and perhaps worthy of further exploration. 

Once this connection is made it is then expected that the the coefficients follow 
the normal law (just as the classical Eulerian numbers do). As a matter of fact, one 
of the proofs (although not the first) of the asymptotic normality of the classical 
Eulerian numbers is via a fairly general device due to Bender ([1]) and nowadays 
referred to as Bender's theorem. This is, indeed, the approach we take and in the 
subsection 6.1.4 we verify the conditions of Bender's theorem to conclude our proof. 
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6.1.1. Bivariate generating function. We give the formula for the joint generating 
function of r„ and A n in the following statement. 

Proposition 10. Let z,t be complex numbers. Then we have 

1 " 

(18) E n t A "z r " = — £(t - l) n - k c n , k (z), 

' fe=0 

where the coefficients {c m j(z) : < £ < m < n} satisfy the following recurrence 

c m +iA z ) = ( z + M)CmA z ) + ( z + 2 ( e ~ !) + l)Cm,«-i(«), 1 < ^ < m, 
wii/i £/ie following boundary conditions: 

c . (z) = 1, c m+ i : „ l+ i(z) = (z + 2m + l)c mjTO (z), £^+1,0(2) = 0C m>o (2:). 

Proof. If n = and we set co,o( 2 ) = 1 then both sides of (18) are 1. Otherwise, let 
Ij indicate the event that we put an a/7 in the jth box on the diagonal (counting 
from NE to SW). Then A n = £)" =1 Ij and we have 

E n t A "z r " = E n E{t A ^ +I -z r -\V n ^) = E n t An - 1 E(t In z Tn |D„_i), 

where we have used the basic properties of the conditional expectations (see [22, 
Formulas (16), p. 79 and (17), p. 80]). Now, I n = 1 means that we put an a/7 in 
the SW corner. In that case we have r n = r„_i + 1 and since this happens with 
probability 1/(2 • 3 r "" 1 ) we get 

E^z^V^) - ^ r "- 1+1 ^ T ^r 7 + E(t ; «/»/;. = o|D B .i). 

The second term is equal to 

E(* r »Jj B=0 |E> n -i) = E(^»|2? n _i)-E(« r »J Jn= i|X> n _i) 
z + l /z + 2\ r "- 1 z r — 1+1 
2 V 3 J ~ 2 • 3 r —i ' 

Combining we obtain 

E ( ^-o^(fr + ^'-'^-' 

Using the second part of (7) leads to the basic recurrence 



= ^ {*(* - l)E„_it yln ~ 1 ^ rn ~ 1 + (z + l)E„_it An_1 (z + 2) r - 1 } . 
Upon further iteration of this relation we obtain for any < m < n: 

1 m 

E n t A -z r - = y"(t-l) m - e c me E n _ m t A "-™{z + 2t) r "-™, 

2 m n(n-l)...(n-m+l) ' V ; 

for some coefhcients c m ^ = c m j(z), < I < m. To see that they satisfy the stated 
recurrence, we apply the basic recurrence (with n — m instead of n and z + 2£ instead 
of z) to the expectation on the right-hand side above. We get that it is equal to 

, 1 J (z + 2£)(t - l)E n _ m _!t A »-™ -Hz + 2^f— 1 
2{n — m) I 

+{z + 2£+ l)E n _ TO _it A "- m - 1 (z + 2(1 + l)) r — 1 }. 
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Substituting this into the above formula for E„i ™z r ™ and multiplying both sides 

by 2 m+1 n(n — 1) (n — m) (to avoid writing a denominator on the right-hand 

side) we get 

m 

J2(t - l) m - e c m ,e{(z + 2£)(t - l)E„_ m _if A "- m - 1 (z + 2£) r — 1 

+(z + 21 + l)E„_ m _if An - m_1 (z + 2(1 + 1 }. 

Splitting this in two sums, isolating the £ = term in the first, and the I = m in 
the second, and then shifting the index in the second sum, we further obtain 

m 

+ 5> - l) m - e c m Az + 2*)(t - l)E„_ m _ii A ' l ~ m ~ 1 (z + 2^) r — 

m— 1 

+ (* - l) m_/ <W(* + 2^ + l)E n _ (m+1) t A »-c«.+D (z + 2(* + 

+c m , m (z + 2m + l)E„_ (ro+1) i A «-<™+i> (z + 2(m + 
= (t - l) m+1 c m , zE n „ (m+1) ^-<™+i)z'-"-('"+i) 

+ J> - l) ro+1 -^c ro ,,(z + 2i)E n _ (m+1) t A -^ (z + 2iyn-im +1) 

m 

+ J^it - l) m " (/_1) cm,«-i(« + 2(* - 1) + l))E„_ (ro+1) i A -(-+^) (z + 2*) r »-<™ 
e=i 

+c mim (z + 2m + l)E n _ {m+1) t A -^ (z + 2(m + . 
So, if we write it as 

m+l 

£ (t - l) m+1 --' Cm+1 ,,E„_ (m+1) ^-<™+ 1 ) (z + 2ff , 

as postulated above, we see that 

Cm+1,0 = ZC m fi, c m +i,m+i = (z + 2m + l)c mjTO 
and, upon combining the two middle sums, that 

c m+ i,e = (z + 2t)c m! z + (z + 2(t-l) + l)c m ,e-i, for \<l<m. 
Taking m = n — 1 we get 

n— 1 

En^ r " - ^ T ^^(t-l)"- 1 - £ C „-i / E 1 ^(^ + 2^r, 

and since 

Eii Al (* + 2*) ri - Ii( z + 2 ^) + ^ = ^(t-l)(z + 2£) + i(z + 2^+l), 
we can write 
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where the coefficients {c m ^ : < t < m < n} satisfy the stated recurrence and the 
boundary conditions. The proof is complete. □ 

6.1.2. Expected value and the variance. A number of properties of (A n ) can be 
deduced from Proposition 10. We illustrate this by obtaining the exact expression 
for the expected value and for the variance of A n . To help facilitate that we put 
the coefficients (c n _k) in a Pascal type triangle. 

co.o 

S \ 

Cl,0 Ci,i 



£n,0 Cn,l ■ • • C-n,n— 1 Cn,n 

The SW move from a coefficient c m _i has weight z + 11 and a SE move has weight 
z + 2£+l. The value of a given coefficient is obtained by summing over all possible 
paths leading to it from the root c 0i o the products of weighs corresponding to the 
moves along the path. For example, there is only one path leading to c n , n (all 
moves are SE) and hence 

n-l 

Cn,n = Y[(z + 2] +1). 

3=0 

Likewise, paths leading to c ni „_i have exactly one SW move; thus there are n of 
them and if the sole SW move is from Ck,k, < k < n — 1, then the weight of that 
path is 

'k-\ \ /n-2 

J](* + 2j + l) (z + 2k) ]J( z + 2j + l) 

^=0 / \j = k 

Consequently, 

n-l /fe-l \ /n-2 

Cn ,„_i(0)=^ n( z+2 ^ +i ) ( z+2fc ) n( z+2 -? +i ) ] • 

k=0 \j=0 J \j=k 

The significance of this is that 

Q 1 
En^n 777 E n £ n Zi — — -C nn _i(l) 

ot t=i,«=i 2™n! 

1 n-l /fc-1 \ /n-2 

' fe=0 \j=0 / \j=k 

= 2 4 ! 2 "-V-l)!E(2 fc + l) = f. 

fe=0 

This proves the first part of (14). Similarly, to compute the variance we use 
var(A„) = EA n (A n - 1) + EA n - {EA n ) 2 and the fact that EA n (A n - 1) = 
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^2 (E£ A ")| t=1 = 2^r c n,™-2(l)- Paths leading to c nj „_ 2 (l) have exactly two SW 
moves, the first could be from any Ck,k, < k < n — 2 and the second from cej-i 
for some k < i < n— 1. These two moves have weights 2fc + l and 2£— 1 respectively, 
and the remaining SE moves have weights 

2 • 1, 2 • 2, . . . , 2k, 2(k + 1), 2(k + 2) . . . , 2{l - 1), 21, . . . , 2{n - 1). 

Therefore, 

n— 2 n— 1 

c„,„- 2 (l) = E E (2fc + l)(2^-l)2"- 2 (n-2)! 

fe=0£=fe+l 



2"- 2 (n - 2)! 2(2* - 1) J2( 2k + ^ 



= 1 fc=0 

-l 



= 2 n - 2 (n-2)!j2(2£- l)^ 2 



Hence 



= 2"- 2 (n-2)!^n(n-l)(3n 2 -5n + l) 

= 2"- 2 n!(in 2 -^n+i). 

v 2 6 6 ; 

, „ , o 2"- 2 n!,l 2 5 1. n n 2 n 2 5 1 6n n 2 n + 1 

var(A>) = 2 ^(-n 2 n+-)-\ = n-\ 1 = . 

V n> 2 n nl V 2 6 6' 2 4 4 12 12 12 4 12 

Thus we have proved the second part of (14) as well. 

6.1.3. Connections to generalized Eulerian numbers. Before establishing the as- 
ymptotic normality of (A n ) we take a closer look at the doubly indexed sequence 
{c n ,k} since it has intriguing connections that may be of interest in their own right. 
It is featured as entry A145901 in [23] and is closely related to another sequence 
from [23], namely A039755. More precisely, c n ^ = 2 k k\W 2 (n,k) where W m (n, k) 
are the Whitney numbers of the second kind satisfying the recurrence 

W m (n, k) = (mk + l)W m (n -l,k) + W m (n — l,k — 1). 

The numbers (W m (n, k)) were introduced in [16] and their properties were studied 
in [2, 3, 4]. Since we are dealing exclusively with the case m = 2 we drop the 
subscript and we write W(n, k) for W2(n, k). Thus, the generating function of A n 
may be written as 

n 

(19) Mt) = £ n t A " = 2 2 fc fcW(n, k)(t - l)"- fc . 

' k=0 

It is perhaps of interest to mention that the numbers (2 k k\W{n, k)) themselves 
satisfy the central (and local) limit theorem as was shown in [5]. However, this is 
not exactly what we want since the generating function above is in powers of t — 1 
rather than t. In terms of powers of t ip n (t) has the following form. 

Proposition 11. The probability generating function of the number o/a/7 entries 
on the diagonal of a random staircase tableau of size n has the form 

1 " 

(20) lM*)=2^iE^ n ' m )* m ' 
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where the numbers {V(n,m), < m < n} satisfy the boundary condition V(n, 0) 
1, £/ie symmetry relation V(n, m) = V(n, n — m), and the recurrence 



(21) V(n, m) = (2m + l)V(n - 1, m) + (2(n - m) + l)V(n - 1, m - 1). 



T/ie explicit expression for V(n, m) is given by 



n—m / , 

'n — k 



(22) V{n,m)=^2 k k\W{n 1 k)\ K )(-l)"-' s - m , < m < n. 



fc=0 



Proof. We start out with (19). To rewrite it in powers of t, note that the mth 
derivative of ip n at t = is 



n—m 

^ m) (°) = ^ E 2 ** 1 ^"' *)(""*) (n-fc-(m-l))(-ir- fe - 



2™n! 

fe=0 



1 n—m / ,n. 

-i- V 2 fc fcW(n,fc) / (n ; kh - (-1) 
2™n! ^ v ; (n-fc-m ! v 7 



n— fc— m 



Therefore, 



m— m— 



which shows that (20) holds with y(n, m) given by (22). 

It remains to verify the claimed properties of the numbers V(n, m), < m < n. 
Since the symmetry condition follows by induction from the recurrence it suffices to 
verify the recurrence and the boundary condition. For the boundary condition, we 
see immediately that for n > V(n,n) = W(n,0) = 1, so that once we verify the 
recurrence (and thus also the symmetry) we will have that V(n, 0) = 1 for all n > 0. 
To verify (21), we use the basic recurrence for W(n, fc)'s to write the left-hand side 
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of (21) as 

k 



V(n,m) = ^2 k k\({2k + l)W{n-l,k) + W{n-l,k-l)M n J(-l) 

fc=0 

n—m / ,\ 

= ^2 k k\(2k + l)W(n-l,k)( n ~ J(-l)"-™- fc 
fe=o \ m J 

n—m / j\ 

+ ^ 2 fc fcW(n-l,fc-l)( n ~ 
fe=i V m / 

n—m / , \ 

= J] 2 fc fc!(2fc + l)VK(n- J(-l) 

n —m— l / _ I. _ 1\ 

+ ^ ^'(fc+Wn-U) " * 1 (-1) 

I — n V / 



n—m—k 



\n—m—k 



k=0 

2 n_m (n - m)!(2(n - m) + l)W(n - 1, n - m) 

n—m—l / , \ 

+ 53 2 fc fc!(2fc + l)T^(n-l,fc)^ m J(-l) 

n—m—l / 

- 53 2 fc fc!(2(fc + l))W(n- l,fc)( 

fc=0 ^ 



n—m—k 



n — k — 1 

m 



(-1) 



n—rn—k 



On the other hand, the right-hand side of (21) is 

n—m—l 



n—m—l / 

(2m + 1) 53 2 k k!W(n-l,k)i 

l — n ^ 



k=0 

n—m 



>; — 1 — k 

m 



(-1) 



n—m—k — 1 



/ b 1 \ 

+(2(n - m) + 1) 53 2 fe fc!VK(n - 1, *) ( " ~ _7 J (-!)"" 



fe=0 

n—m—l 

\n—m—k — l 



(2m + 1) 53 2 k k\W{n- l,fc)(™ & 1 )(-1) 



m—k 



k=0 

n—m — l , 7 -i \ 

+(2(n-m) + l) 53 2 fc fcW(n-l,/fc)( n ^_~ V-l)"-"*-* 
fe=o \ m / 

+(2(n - m) + l)2"- m (n - ro)!W(n - 1, n - m). 

So, we see that the coefficients in front of W(n— 1, n—m) in both expressions are the 
same, and to complete the verification of (21) we need to see that the coefficients 
are the same for the remaining values of k, < k < n — m — l. Cancelling the 
common factors, we need to see that 

,-,fn—k\ ^.fn — k — 1 

(2fc + l) -2(fc + l) 

\ m J \ m 

,[n — k — l\ . , / n — fc — 1 

= (2(n-m + 1 1 -(2m + l) 

\ m — 1 / V m 
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But that is straightforward: using (™ m fe ) = ( n 1 ) + (™ m 'l 1 1 ) and grouping the 
terms this boils down to verifying that 

" ~ k ~ 1 )(2k + 1 - 2(k + 1) + 2m + 1) = ( U ~ k ~ (2(n - m) + 1 - 2k - 1), 
m J \ m — 1 / 

or, equivalcntly, that 

/n — fc— 1\ , ,/n— fc— 1 

to = (n — to — fc) 

\ m J \ to — 1 

which follows immediately from the defining property of the binomial coefficients. 

□ 

Remark: Triangle of numbers F(n, to), < to < n is featured in the Online 
Encyclopedia of Integer Sequences [23] as a sequence A060187 (with the shift in 
indexing: V(n, to) = T(n + 1, m + 1)) and is called "Eulerian numbers of type B" . 
This sequence can be traced back in the literature to MacMahon's paper [21] and 
it was subsequently studied in more detail in [19, Sec. 3.2] as a sequence B n ^(l). 
In particular, it appears that the expression for a bivariate generating function of 
(V(n, k)) was derived for the first time in [19]. We use this expression in the next 
section to derive the asymptotic normality of (A n ). 

6.1.4. Conclusion of the proof by Bender's theorem. Using the properties of the 
numbers V(n, to) given in Proposition 11 (and the identification with the sequence 
A060187 from [23]) we can complete the proof of Theorem 9 by establishing (15). 
To do that we will rely on a general theorem due to Bender [1, Theorem 1]. (Ben- 
der result is also described in [18]; see Section IX. 6 in general, and Theorem IX. 9, 
Example IX. 12, and Proposition IX. 9 in particular). Recall from [23] or [19, For- 
mula (3.23)] (and see Section 3 of [19] for a proof) that the bivariate generating 
function of the numbers (V(n, k)), called in [19] (B ni k(l)), is 

w k z n (1 - w)e^' w >- 



n>0 k=0 

Therefore, it follows from (20) that the bivariate probability generating function of 
the sequence (A n ) is 

w k /z\ n (1 — i/iVC 1 -™) 2 / 2 

n>0 n>0fc=0 

where we define f(z, 1) = 1/(1 — z). We now closely follow the way Bender applied 
his result. First, 

J 6 > l _ e s e (l-e")z 

has a simple pole at z = r(s) = s/(e s — 1). Furthermore, 

e (l-e s )2 = e (l-e 3 )( Z -r( S )) e (l-e=)r( S ) = e (l-e 3 ) (z-r(s)) g - S ^ 

so that 

_ (1 _ e «) e -«/2 e (i-e')(*-r( fl ))/2 _ (1 - e s )e- s ' 2 

^ Z,e ' ~ 1 _ e (l-e»)(*-r(«)) - e ( e .-l)(z-r( S ))/2 _ e -(e'-l)(z-r( fl ))/2 ' 

Since for bounded u _ 

e u - e- u ~ 2u + 0<y1 ^ 
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we get 

fl - e s )e~ s/2 
(e s — l)(z — r(s)) 

with the constant in 0(1) bounded when both s and z — r(s) are close to 0. Thus, 
by a comment at the very beginning of Section 3 of [1], the conditions of Theo- 
rem 1 of that paper are satisfied, and hence the central limit theorem holds (with 
centering by EA n ~ n/2 and scaling by y/vax(A n ) ~ y/n/12 as implied by (14)). 
Alternatively, we see that 

r(0) = l, r'(0) = -i, and r"(0) = ^ 

so that by [1, Theorem 1] EA n ~ ^ and var(A„) <~ j| which conforms to what we 
have found in (14). In any event, (15) follows. 



7. Conclusion and further remarks 

In this paper we have developed a probabilistic approach to the analysis of 
properties of random staircase tableaux. Using this approach we established the 
asymptotic normality of several parameters associated with appearances of Greek 
letters a, (3, 7, and S in a randomly chosen tableau. We certainly hope that this 
approach will be useful in the analysis of other properties of staircase tableaux. 
From the combinatorial point of view, it would be of interest to analyse the number 
of appearances of the letter q in such tableaux. It is not clear at this point, that 
the method we develop is adequate to address that question, and if so, how difficult 
it would be to achieve. This is probably an issue worth resolving in the future. 
The fact that our method can be used to give new and rather complete results 
concerning Greek letters and the fact that it gives an easy way of enumerating of 
staircase tableaux of a given size makes us cautiously optimistic. 

It is perhaps worth making the following point. In our arguments we relied 
on observations like (17) and symmetries of staircase tableaux to minimize the 
amount of work. We wish to emphasize, however, that our probabilistic approach 
does provide a unified and systematic way of analyzing each of the statistics we 
considered in a self-contained manner (i.e. not relying on relations between various 
parameters). In fact, a direct proof that T n satisfies (11), (12), and (13) was given 
in [13]. This is partly a reason we believe that our approach has potential of being 
useful in addressing other questions concerning statistics of staircase tableaux. To 
reiterate that point we wish to briefly sketch a proof of (16) in Theorem 9 not 
relying on symmetries of the parameters. 



7.1. The number of 0/6 on the diagonal. We show that A n and B n have 
identical probability generating functions which, of course, implies (16). Write 
B n = Y^j=i ty, where bj is 1 if we put a j3 or 5 in the jth box on the diagonal 
(counting SW from the top) and is otherwise. We first derive the expression for 
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the conditional probability generating function: For z, t complex, 

E(t 6 »z r »|2? n _i) - £{t b -z r -I bn=l \D n ^) + E(t h -z r -I K ^\D n ^) 
= tE(z r -I bn=1 \V n ^) + z r ^ +1 P(b n = 0|X> n _i) 

= iE(* p »|2>„_i) - fE(z r "/ 6 „ =0 |X>„-i) + « r »- 1+1 



= f 



' +i (^r-^-^^(sF ' 



2 



This gives an expression for the bivariate probability generating function for r„ and 
£?„ and a recurrence for the coefficients: 

^ = ^-.{ t i±i(i±2)-' + «,-«,|(§r} 

= ^ {*(* + lJEn-it^"" 1 (z + 2)^ + (1 - t)zE n _ 1 t B «- 1 z r »- 1 } 

= ^7 ^7 —\y j b m ,eE n _ m t B ^{z + 2i)^-A. 

2n2(n — 1) 2(n — m+1) 1^ j 

Now, 

E n _ rn t B »-»(« + 2^) r »-» = — ,{i(z + 2£ + l)E„_ m _ 1 t s — !(z + 2^+l) r — 1 

+(1 -*)(« + 2^)E„_ m _ 1 t B — 1 (z + 2^) r "-— 1 }, 

so that 

6 mi£ E„_ m i B — (z + 2£) r — = t(« + 2£ + l)6 m ,,E n _ m _ 1 t B '— 1 (z + 2(£ + 1))^— i 

+(1 -t)(z + 2£)b m , i E n _ m _ 1 t B —- 1 (z + , 
which means that the coefficient in front of 

E n _ (m+1) t B -^(z + 2£)^-(^ 

is 

(1 -t)(z + 2£)b m , e + t(z + 2{l - 1) + l)6 m ,,_i. 

So, if we write b m ^ = a m jt l {l — i) m ~ l we see that the coefficients a m ^ — a m j{z) 
satisfy the recurrence 

fflm+M = ( z + 2 t) a m,i + (z + 2£- l)a m> e-i, 

with the initial condition 

ao.o = 1, and a m j = 0, for £ > m. 

This is exactly the same recurrence as for the sequence (c n ^) defined in the last 
section and thus 

n 

(23) E n t B *z r « = —J2c n , k t k (l-t) n - k . 

' k=0 
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In particular, putting z = 1, differentiating with respect to t, evaluating at t = 1, 
and using c„^„ = 2™n! and c„ jn _i = 2 n ~ 1 n\n we confirm that 

1 1 n 

EnB n = — ,(-c n ,„_i +nc n ,„) = — — (-2"- 1 n!n + n2"n!) = -. 
2"n! 2™n! 2 

Finally to see that B n and ^4„ have, in fact, the same distribution, put z = 1 in 

(23) and expand (1 - t)"- fe to get 

i-k 



5>„, fei fe (i-t)"- fe = • (- 1 )" Hr '" i 

fc=0 fe=0 j=0 \ / 

n n— k / , \ 

fc=0j=0 \ J / 

= E'™{e^(;::)(-i)™-*} 

m=0 U=0 V 7 J 



m=0 U=0 



Now recall that by symmetry 



rt-(n-m) , _ t\ 

V(n,m) = V(n,n-m)= ^ 2 fc fcW(n, k) I J (_ 1 )n-fc-(n-r„, 



fc=0 

so that 



m— fc 



fe=0 m=0 

which is exactly the same as the generating function of A n as claimed. 
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